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ABSTRACT 

We discuss the structural and morphological properties of galaxies in a z = 1 .62 proto-cluster using near-IR 
imaging data from Hubble Space Telescope Wide Field Camera 3 data of the Cosmic Assembly Near-IR Deep 
Extragalactic Legacy Survey (CANDELS). The cluster galaxies exhibit a clear color-morphology relation: 
galaxies with colors of quiescent stellar populations generally have morphologies consistent with spheroids, 
and galaxies with colors consistent with ongoing star formation have disk-like and irregular morphologies. 
The size distribution of the quiescent cluster galaxies shows a deficit of compact (< 1 kpc), massive galaxies 
compared to CANDELS field galaxies at z = 1 .6. As a result the cluster quiescent galaxies have larger average 
effective sizes compared to field galaxies at fixed mass at greater than 90% significance. Combined with 
data from the literature, the size evolution of quiescent cluster galaxies is relatively slow from z — 1 .6 to 
the present, growing as (1 +z)~°' 6±01 . If this result is generalizable, then it implies that physical processes 
associated with the denser cluster region seems to have caused accelerated size growth in quiescent galaxies 
prior to z = 1.6 and slower subsequent growth at z < 1.6 compared to galaxies in the lower density field. The 
quiescent cluster galaxies at z = 1 -6 have higher ellipticities compared to lower redshift samples at fixed mass, 
and their surface-brightness profiles suggest that they contain extended stellar disks. We argue the cluster 
galaxies require dissipationless (i.e., gas-poor or "dry") mergers to reorganize the disk material and to match 
the relations for ellipticity, stellar mass, size, and color of early-type galaxies in z < 1 clusters. 
Subject headings: galaxies: clusters: general — galaxies: clusters: individual (XMM-LSS 02182-05102) — 
galaxies: evolution — galaxies: elliptical and lenticular, cD — galaxies: high-redshift — 
galaxies: structure 



1. INTRODUCTION 

Massive elliptical and early-type galaxies dominate re- 
gions of high density suc h as those of galaxy clusters in 
the pr esent Universe (e.g.. lDresslerlll980l iPostman & Gellerl 
119841) . By z < 1.5, Hubble Space Telescope (HST) observa- 
tions show that these passive cluster galaxies have elliptical 
and lenticular morphologies, with a strong color-density re- 
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lation ship (e.g.. Ivan pokkum & Franxll2001 
2001 IPostman et al.1 1200 51: iMei et al.l l2006 
2006tlHilton et al.ll2009t) . The emerging picture for formation 
and evolution of the massive, red, early-type cluster galax- 
ies is one in which these galaxies forme d their stars at z > 2 , 
with subsequent passive evolution (e.g., St anford et al.ll 1998k 



lEisenhardt etal] l200l IWhilev et al.l l200l) . These galaxies 
continue to grow by mergers and secular processes, with neg- 
ligible additional star-formation in order for their color evolu- 
tion to be consistent with observations. 

The details of this evolution is unclear, yet these cluster 
galaxies must assemble sometime. It may be that the for- 
mation of cluster gala xies is related t o the cluster assem- 
bly process itself (e.g.. lDubins ki 1998). Observations show 
that the intracluster galaxy velocity dispersion is lower in 
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mergers, which lead to a population of spherical, elliptical 
galaxies dNavarrol 119901) . Out to z ~ 1, ellipt icals domi- 
nate the galaxy populations of massive clust ers (Des ai etafl 
I2007t iHolden et al.ll2009t IVulcani et al.ll201 ll) . while lenticu- 
lar and early-type spir al galaxies dominate the cores of som e 
lower density groups (IWilman et al.ll200"9t I Just et al.ll2.010h . 
If lower mass groups are common precursors to galaxy clus- 
ters, then their galaxies must undergo morphological evo- 
lution to early-type galaxies as the groups merge to form 
larger clusters. This is expected based on some semiana- 
lytic models, which predict that processes associated with the 
cluster fo rmation are expected to influence galaxy evolution 
at z » 1 dDubinskilll998h iLin & Moh ril2004t iDe Lucia et alj 
l2007tlRuszkowski & Springelll2009l) where the main progen - 



itors of clusters collapse (e.g., iBovlan-Kolchin et al.1 l2009h . 
Therefore, if this hypothesis is correct, then as we encroach 
on the formation epochs of today's massive clusters, z > 1.5, 
we should expect to see rapid evolution in the properties of 
the cluster galaxies. 

In addition, observations show quiescent galaxies (not only 
those in clusters), with apparent early-type morphologies, 
undergo strong size evolution with redshift out to z = 2 (e.g. , 
Daddi et al.l 12005b IPapovich et al.ll200l: iTruiillo et "all 120061 
2007tlLonghetti et all2007tlZirm et al.l2007l:lToft et al.ll2007l: 
Buitrago et all l2008t ICimatti et al.l l2008t Ivan der Wei et al.l 
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Cassata et alj 120101: ISaracco et al.l 120101) . One explanation 



for this size evolution is that these galaxies grow by frequent 
dissip ationless (i.e., gas-p o or or "dry") mino r mergers 
(e.g. lLoeb & Peebles! 120031: Ivan Dokkuml 12005b [Bell et all 
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Lotz et al.1 120081: iMasiedi et alj 120081: i Hopkins et al.l l2009t , 
20 10t Ivan der Wei et al.ll2009al 1201 ll) Minor mergers would 



cause the galaxies to add mass at larger radii, increasing 
their effective sizes su bstantially with a relatively small 
increase in stellar mass (lOser et al.ll2010l). Some recent ob- 
servations support this interpretation ( Bezanson et al.l l2009t 
IHopkins et ai1 l2009at Ivan Dokkum et al.ll2010l) . although this 
explanation would not explain the substantially larger central 
densities of high-redshi ft ellipticals compare d to galaxies at 
lower re d shift ( see e.g., IStockton et al.ll20Tol) . Alternatively, 
iGrahaml (1201 ll) notes that many of these compact objects 
share sizes, masses, and mass densities of present-day bulges, 
suggesting some of these objects are the precursors to the 
spheroidal components of present-day disk galaxies. 

It is unclear how the assembly of ellipticals in high den- 
sity cluster (and forming cluster) regions differs from that 
in the lower density the field. If the size growth of ellip- 
ticals is driven by minor mergers and galaxies experience 
more mergers in forming clusters, then it follows that the 
size and morphological evolution of cluster ellipti cals should 
be acc elerated during the cluster formation stage. IStott et all 
(1201 ll) report that the sizes of the most massive galaxies in 
clusters increase by a t most 30% during the period of z = 1 
to 0.2. ICooper et alj d201 ll) find a correlation between the 
sizes and local galaxy overdensity for early-type galaxies at 
0.4 < z < 1.2, suggesting accelera ted morpho l ogical evolu- 
tion in regions in higher density. IZirm et aTl (1201 ll) find a 
hint of evidence that massive quiescent galaxies in the vicin- 
ity radio galaxy MRC 1138-262 at z = 2.2 have larger sizes 
at fixed mass compared to galaxies in the field at this red- 
shift. These observations support the hypothesis that cluster 
ellipticals experience accelerated structural evolution. How- 
ever, other observations at higher redshift (z ~ 2.3-4.1) find 



no evidence that the sizes or morphologies of galaxies dif- 
fer in high- density regions compared to those of lo w-density 
regions (e.g. iPeter et alj|2007t lOverzier et aT]|2008l) . suggest- 
ing that any environmental effects are not yet present at these 
epochs. 

Here, we compare the properties of galaxies in the high 
density region of a f orming cluster at z = 1 .62, XMM- 
LSS J02182-05102 (IPapovich et al.1 l2010t iTanaka et alj 
l2O"l0l) . and we compare them to similarly selected galaxies in 
the lower density z = 1 .6 field. This galaxy cluster was identi- 
fied as an overdensity of sources with Spitzer /IRAC colors in- 
dicative of high-redshift galaxies (IPapov ich 2008). The clus- 
ter shows a dominant population of red galaxies, which form a 
strong "red sequence" population, with an esti mate of the last 
major star-formation epoch ofz/ = 2.2-2.3 ( IPapovich et alj 
120101) . In addition, this cluster shows a significant fraction 
of star-fo rming galaxies as evidenced by their SpitzerllA /im 
emission dTran et al.l 120101) . There are currently 13 redshifts 
for galaxies with 1.62 < z < 1.65 within a physical projected 
radius on the sky of 1 Mpc of the cluster center (10 of these 
galaxies have 1 .62 < z < 1 .63; Papovich et al. 2010, Tanaka et 
al. 2010, 1. Momcheva et al., in prep, C. N. A. Willmer et al., 
in prep). These redshifts provide an estimate of the velocity 
dispersion and total cluster mass assuming the cluster is viri- 
alized, ay = 360 km s" 1 and M 2 oo ~2x 10 13 M Q , although 
there is evide nce to suggest the assumption of yirializ ation is 
unlikely (see IPapovich et af]l20 1 Ot iPierre et al JI20 1 ll) . which 
is entirely consistent with the expected assembly histories 
of a present-day massive c luster observed at z ~ 1.6 (e.g., 
IBovlan-Kolchin et al.ll2009l) . While the reported velocity dis- 
persion was consistent with the weak (4a) XMM X-ray de- 
tection (IPapovich et al.|[20Tol) . recent Chandra data show that 
several point sources dominate the X-ray emission with very 
faint extended emission, supporti ng the interpretati on that this 
cluster is in the act of collapsing (IPierre et al.ll201 ll) . 

In terms of semantics, throughout this paper we refer to 
XMM-LSS J02182-05102 as a "cluster" even though it is un- 
likely to fully satisfy the classical definition of a virialized 
object. The distinction "proto-cluster" or "forming cluster" is 
strictly more apt as it seems likely that this structure is in the 
process of collapse and assembly. Regardless, because this 
object corresponds to a clear high surface density o f galaxies 
at z = 1 .62 (20(7 as defined by IPapovich et alJl2010l) . we have 
the ability to compare and contrast the morphological evolu- 
tion of galaxies in a high density region compared to that in 
the lower density field. 

The outline for this paper is the following. In § 2 we de- 
scribe the properties of the imaging datasets and we describe 
our analysis. In § 3 we discuss the color-morphology re- 
lation in this cluster. In § 4, we discuss the size-mass re- 
lation for quiescent galaxies associated with the cluster and 
compare it to a similarly-selected sample in the field. In 
§ 5 we discuss the distributions of ellipticities and surface- 
brightness profiles for the quiescent galaxies in both the clus- 
ter and field. In § 6 we consider possible evolutionary sce- 
narios for the quiescent galaxy population, and we discuss 
how environmental processes affect the galaxies' evolution. 
In § 7 we summarize our conclusions. Throughout this pa- 
per we report magnitud es measured relative to the AB sys- 
tem (lOke & Gunnlll 983*1) . We denote photometric magnitudes 
measured in the WFC3 F125W and F160W passbands as /125 
and //i6o, respectively. Throughout we assume a cosmology 
with Q m = 0.3, n A = 0.7, and H = 70 km s _1 Mpc -1 . 
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2. DATA AND ANALYSIS 

The z = 1 .62 cluster XMM-LSS J02 1 82-05 102 is located in 
the UKIRT IR Deep Sky Survey (UKIDSS: ILawrence et all 
120071) Ultradeep survey (UDS). This cluster received par- 
tial HST/WFC3 imaging in the F110W and F160W bands 
as part of the Cosmic Assembly Near-IR Deep Extragalac- 
tic Legacy Survey (CANDELS) program (Pis: S. Faber, 

H. FergusonJJ The CANDELS strat egy, data acquisition , and 
data reduction are descr ibed fully in lGrogin et al.l (120111) and 
iKoekemoer etail j20Tlb . The CANDELS imaging achieves 
limiting magnitudes of J125 = H^o = 26.6 mag (10er for aper- 
tures of 0."4 diameter). Owing to the CANDELS field place- 
ment, the H5T/WFC3 imaging covers slightly more than 50% 
of the galaxies associated with the z = 1 .62 cluster (see below), 
including 6 galaxies with spectroscopic redshifts 1 .62 < z < 

I. 65 within a physical projected radius of 1 Mpc of the clus- 
ter center. The CANDELS imaging does cover most of the 
cluster core including its most massive, quiescent galaxies. 

In addition to the HST imaging, this field has deep 
BRiz imaging f r om th e Subaru-XMM Deep Survey (SXDF; 
[ Furusawa etail 120081) . JK imaging from UKIDSS (e.g., 
iWilliams et al.ll2009l) . Spitzer IRAC data in four bands prob- 
ing 3.6 to 8.0 pm, and MIPS data at 24 /iiriQ 

2.1. Merged Catalogs, Photometric Redshifts, and Sample 
Selection 

As in lPapovich et al.l (120101) . w e used the AT-band sele cted, 
SXDF and UDS catalogs from IWilliams et ail (120091) and 
merged these with the Spitzer/IRAC data. Following Pa- 
povich et al., we used the muliwavelength photometry to de- 
rive photometric redshift probability distribution functions , 
P(z), for each source using EAZY dBrammer et al.l [2008h . 
Here we considered a sample with J < 24.5 mag, which 
is an approximate 3 a limit for the UKIDSS data. As in 
IPapovich et ail(l20Tol) . we define a likelihood that galaxies are 
associated with the cluster redshift, 

V z = J P(z)dz, (1) 

integrated over the redshift range given by z = z cen ± Sz with 
Zcen = 1 -625 and 8z = 0.05 x (1 + z cen ), approximately the 68% 
confidence range on the photometric redshifts for the red, qui- 
escent galaxies. 

We consider all galaxies with V z > 0.3 and projected dis- 
tances R pro j < 1.5 Mpc to be associated with the cluster. 
Galaxies with well-established spectral features, such as the 
4000A/Balmer break, have sharp P{z) and thus higher V v 
which includes red galaxies with lower implied specific star- 
formation rates. Galaxies that are actively star-forming have 
weaker 4000A/Balmer breaks, have more broad P{z), and 
have lower V z . Therefore, choosing V z > 0.3 ensures that 
we do not bias ourselves away from the (bluer ) star-forming 
objects (see discussion in lPapovich et al.ll"2010l) . 

However, in the sections § 4 and 5, we focus on the prop- 
erties of a sample of quiescent galaxies in the cluster to those 
in the field. For this sample of quiescent galaxies we increase 
our selection criterion to V z > 0.5. We do this because the 
quiescent galaxies have tighter P(z) functions, and will have 
higher V z . Our tests have shown that a V z > 0.5 criterion 
provides a cleaner sample as the samples would otherwise 

1 http://candels.ucolick.org/ 

2 http://irsa.ipac.caltech.edu/data/SPITZER/SpUDS 



include galaxies with more than 50% of their P(z) outside 
the desired redshift range. From this subsample we define 
quiescent galaxies associated with the cluster as those with 
^proj < 1.5 Mpc, and those in the field as R pro j > 3.0 Mpc. 

Table [TJlists the properties of the objects in the CANDELS 
cluster and field samples including the astrometric coordi- 
nates, magnitudes, colors, and photometric redshift informa- 
tion. The table includes all galaxies at z = 1.6 satisfying 
V z > 0.3 as defined above. 

2.2. Stellar Masses 

We fitted the 10-band galaxy photometry covering 0.4- 
8 pm with model spectral energy distributions to estimate 
the stellar masses for the galax ies in the sample using the 
method of lPapovich et al.l ( 120011) . We used models for a range 
of stel lar population properties from the iBruzual & Chariot! 
( 120031) stellar population sy nthesis model s allow ing for a 
range of extinction using the ICalzetti e~a l. (2000) law. We 
opt to use the 2003 version of the Bruzual & Chariot mod- 
els to faci l itate t he comparison to other studies, including 
iShen et aTl (120031) . Our tests showed that using the 2007 up- 
dated version of the Bruzual & Chariot models yields stel- 
lar masses systematically lower by a 0.2-0.3 dex, but this 
does not affect our conclusions. We assumed models with 
solar metallicity and a Chabrier initial mass function (IMF; 
using a Salpeter IMF would to first order increase system- 
atically the stellar masses by ~0.27 dex). Given that most 
of the galaxies a ssociated with the cluster are quite massive 
dTran et aHl2010[). the solar-metallicity a ssumption is reason- 
able (see, e.g.. Ivan Dokkum et al.ll2004l) . Using different as- 
sumptions for the stellar population m etallicities will affect 
the derived stellar mas ses by ~0.2 dex ( IPapovich et al.ll200Tt 
iMarchesmT et al. 2009). We generate a multi-parameter prob- 
ability distribution function for each galaxy from this mod- 
eling. We then compute the mean and 68% confidence re- 
gion on the stellar mass for each galaxy by ma rginalizing over 
the ot her model parameters (see discussion in IPapovich et al.l 
120061) . Our analysis of the spectral energy distributions of 
galaxies provides an estimate of the instantaneous SFR, which 
we measure as the SFR averaged over the prior 100 Myr using 
the best-fit stellar-population model. Table[T]lists the derived 
stellar masses and SFRs for each object in the sample. 

2.3. Galaxy Morphologies and Sizes 

We used GALFIT dPeng et al.ll2002b to fit models to galax- 
ies in the CANDELS WFC3 F125W imaging, from which we 
determined effective radii, R e ff, and Sersic indices, n. The 
models assume the surface brightness of the galaxies is pro- 
portional to exp(-R/R ef f) l / n dSersidl 1968b . where R is the an- 
gular radius from the galaxy center and where the Sersic index 
is a c oncentration parameter . An exponential disk has n = 1, 
and a ide Vaucouleursld 19481) profile has n = 4. GALFIT con- 
volves the models by the image PSF before fitting them to the 
data. We generated model PSFs for each dither position and 
orien tation for th e two WFC3 imaging epochs using Tiny Tim 
v7.2 (lKristlll9"95l) . The PSF models were dithered and com- 
bined in the same way as the CANDELS data. 

We fitted each galaxy with GALFIT, keeping the position, 
background, orientation, effective semimajor axis, Sersic in- 
dex, and ellipticity as free parameters. We used the WFC3 
F125W image for this analysis as this bandpass corresponds 
approximately to the rest-frame B-band at z = 1 -6, facilitating 
the comparison to other datasets. Our tests show that none of 
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J - [3.6] (mag) 

Figure 1. The observed J — [3.6] versus z— J color-color diagram for galax- 
ies at z = 1.6 in the CANDELS UDS field. At z = 1.6 these observed colors 
correspond approximately to V-J and U — V rest-fram e, which is very effec - 
tive in separating quiescent and star-forming galaxies I Williams et aU2009t) . 
Quiescent galaxies are expected to populate the upper left region of the plot 
denoted by the polygon. Star-forming galaxies form a sequence below the 
quiescent region, where the arrow illustrates the expected change in color for 
A(V) = 1 mag of dust extinction for a galaxy at z = 1.6. The symbol colors 
scale with the specific SFR (sSFR, the SFR per unit stellar mass) as defined 
by the inset color bar. Small yellow circles denote 24 /jm-detected sources 
with /„(24/zm) > 50 /ijy. 

our conclusions would be strongly affected if we instead used 
the WFC3 F160W image. During the analysis we required 
n < 6 because higher values of n usually do not improve the 
fit, and the covariance between n and the effective radius leads 
to an overestimate of R e ff for larger «. Only three objects in 
our samples had best fits with n > 6, and we refit those objects 
forcing n = 4. In what follows, we analyze the galaxies' ellip- 
ticities, defined as e = 1 -q, where q = b/a is the ratio of the 
semiminor to semimajor axes calculated by GALFIT. Table Q] 
lists the GALFIT measurements for all objects in the samples. 

The effective sizes we report in this paper are the circu- 
larized effective radii, R e g = \fab~ = a e f[y/q, where a e ff is the 
effective semimajor axis measured by GALFIT and other val- 
ues are as above. The circularized effective radius is smaller 
than the effective semimajor axis, but it is commonly used 
in the literature, and we use it here for comparison. The 
effective semimajor axes can be computed using the infor- 
mation in Table |T| Furthermore, we have checked that the 
circularized effective radii from GALFIT are in good agree- 
ment with independent, non-parametric measurements of the 
galaxy half-ligh t radii computed following the methods in 
lLotzetalJd2008h . 

We performed a series of simulations to estimate the errors 
in the GALFIT parameters. We inserted model galaxies of 
known effective radius, Sersic index, and magnitude into the 
WFC3 F125W data, and we recovered their parameters us- 
ing GALFIT as desc ribed above. As with other studies (e.g., 
lHaussler et al J 12007). we find that the errors in effective ra- 
dius and Sersic index are correlated strongly, with larger un- 
certainties on effective radius for objects with larger Sersic 
indices. Quantitatively, our simulations show that the mea- 
sured effective radii are accurate to better than 40% for sim- 



ulated compact objects (measured n = 4 and r e ff < 0.5 arcsec) 
with magnitudes typical of the faintest objects in our samples, 
ra(F125W)=23 mag. Similarly, the measured Sersic indices 
are accurate to better than 20% for these objects. The un- 
certainties are substantially lower for brighter and less com- 
pact ( n < 4) objects (similar to the findings of lHaussler et alJ 
120071) . In practice, the errors on the Sersic index have no sub- 
stantive impact on our conclusions. 

3. THE COLOR-MORPHOLOGY RELATION IN A Z = 1.62 CLUSTER 

We select quiescent galaxies using a J -[3.6] and z-J 
color-color selection. At z = 1.6 these colors correspond ap- 
proximately to rest-fr ame V-J versus U - V at z = 1.6, which 
IWilliams et al.l (120091) showed effectively separa tes quiescent 
galax ies from star-for ming galaxies (see also IWuyts et"aT] 
12009b iPatel eTaT]l201 ll lOuadri et al.ll201 ll) . Figure ffl shows 
the J- [3.6] versus z-J diagram for all galaxies in the CAN- 
DELS UDS field with with V z > 0.3. The symbol colors de- 
note the specific SFR (sSFR: the SFR per unit stellar mass). 
Star-forming galaxies form a sequence below the quiescent 
galaxies where the slope of the sequence roughly follows the 
expected change in color associated with dust extinction. 

Quiescent galaxies lie in the upper left region of figure[H as 
indicated by the polygon defined by 

(z-J)ab> 1.3 mag 
(7-[3.6]) AB <2.1mag (2) 
(z-J)ab> 0.5 + 0.55(7- [3. 6])ab. 

We define samples of quiescent galaxies as satisfying all the 
color criteria of equation [2] Based on the analysis of the 
galaxies' spectral energy distributions (§ 2.2), galaxies se- 
lected using the color selection above in the CANDELS sam- 
ple have low specific SFRs. We find that 69/78 (88%) of the 
quiescent galaxies have specific SFRs < 10~ 2 Gyr (includ- 
ing all but one of the quiescent galaxies associated with the 
Z = 1 .62 cluster). Therefore, the quiescent galaxies selected by 
the color selection above have spectral energy distributio ns in- 
dicative of highly "suppressed" SFRs ( iKriek et al.ll2006l) . The 
MIPS 24 /im data give an independent measure of star forma- 
tion or the presence of an AGN. Few of the quiescent galaxies 
are detected at 24 /im: only 3 out of 24 cluster galaxies and 
6 out of 72 field galaxies have /„(24yum) > 50 fiJy. We do 
not reject these sources from the quiescent sample because 
the source of the 24 ^m emission in these galaxies is uncer- 
tain. However, given the small number of 24 /im sources, our 
tests show that none of our conclusions would change if we 
did remove these sources. 

The galaxies associated with the z = 1 .62 cluster exhibit 
a clear color-morphology relation. Figure [2] shows the HST 
WFC3 (F125 and F160W) color images for the cluster galax- 
ies with V z > 0.3 and projected distances R < 1.5 Mpc from 
the cluster center. Spheroids dominate the morphologies of 
the cluster galaxies with colors of quiescent stellar popula- 
tions defined by equation [2] In most cases these galaxies are 
highly symmetric with elliptical and lenticular morphologies 
and a range of sizes. Galaxies with colors consistent with on- 
going star formation have disk-like and irregular morpholo- 
gies. In several cases galaxies in the star-forming region show 
evidence for multiple components, including apparent bulge 
and disk morphologies. This is especially visible in the star- 
forming galaxies with redder/- [3.6] colors, and these galax- 
ies appear to have large effective sizes compared to the bluer 
star-forming galaxies. 
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Figure 2. The observed J — [3.6] versus z—J color-color diagram for galaxies associated with the z = 1.62 cluster. The plot includes all objects within 1.5 Mpc 
(projected) of the cluster with V- > 0.3 (see text) and that have HST coverage from CANDELS. The color images show 6" X 6" (approximately 50 kpc X 
50 kpc at z = 1.6) cutouts from the CANDELS WFC3 F125W and F160W data. The images are placed at the approximate measured color each galaxy (slight 
adjustments to the measured colors have been applied for clarity, but these shifts have no affect on the conclusions). There is a clear color-morphology relation 
in the galaxies associated with this cluster. 



Figure [3] shows the same J - [3.6] versus z—J color-color 
plot as in figure|2]with the galaxies denoted by symbols based 
on their Sersic indice s as measured by GALFIT. Motivated 
by iHogg et all (120041) . we classify galaxies with high Sersic 
index, n > 2, low Sersic index, 0.5 < n < 2, and very low 
Sersic index n < 0.5. Galaxies show a relation between their 
Sersic indices and their location in the color-color plot of fig- 
ure [3] Most of the galaxies with high Sersic indices fall in 
the region of the plot occupied by quiescent galaxies: the qui- 
escent galaxies have surface-brightness profiles dominated by 
spheroids. Galaxies with low and very low Sersic indices fall 
primarily in the region of the plot occupied by star-forming 
galaxies: they have surface-brightness profiles dominated by 
disks. Quantifying these statements, we find that of the 24 



cluster galaxies with z — J and J- [3.6] colors of quiescent 
galaxies, 19 (79%) have n > 2, suggesting a high early-type 
galaxy fraction among the passive galaxies in the cluster. Of 
the galaxies in the star-forming region of figure [3] 28 of 38 
galaxies (74%) have n < 2, implying they are dominated by 
objects with disk-like or irregular morphologies. Further- 
more, based on the simulations in § 2.3 the errors on the Sersic 
index have no substantive impact on our conclusions. There- 
fore, the color-morphology relation exists in this z= 1.62 
cluster, with high Sersicindex (spheroid-dominated) galaxies 
populating the quiescent region of the color-color plot, and 
with low Sersicindex galaxies populating the star-forming re- 
gion. This extends a si milar result obser ved for field galaxies 
dWuvts etal1l2011l) and dBell et alj|2011l) to higher density re- 
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Figure 3. J— [3.6] versus z—J color-color plot for galaxies associated with 
the z — 1.6 cluster. The symbols denote the galaxy Sersic indices, n, as la- 
beled. As in rigure|2] the plot includes all objects within 1.5 Mpc (projected) 
of the cluster with V z > 0.3 (see text) and that have HST coverage in CAN- 
DELS. Also as in figure [2] the data points are placed at the approximate 
measured color each galaxy. There is a clear relation between the galax- 
ies' morphological Sersic index and their location in the color-color plane. 
The galaxies in the quiescent region of the plot have high Sersic indices, in- 
dicative of galaxies with spheroid-dominated morphologies. Galaxies with 
colors of star-forming galaxies have lower Sersic indices, indicative of disks 
and irregulars. Sources denoted by gray diamonds have 24 fim detections 
with f„(24fim) > 50 /ijy. 

gions associated with the cluster at these redshifts. 

4. THE SIZE-MASS RELATION FOR QUIESCENT GALAXIES AT 

Z = 1.6 

4.1. Comparison between Cluster and Field Quiescent 
Galaxies 

Figure |4] shows the (circularized) effective radii of the the 
quiescent galaxy samples from CANDELS in both the z = 
1.62 cluster and the z = 1.6 field as a function of their stel- 
lar mass. At fixed stellar mass, the quiescent field galax- 
ies at z = 1 .6 in the CANDELS data generally have effec- 
tive radii smaller by about a factor of 3 compared to the 
distribution of low redshift early- type galaxies from SDSS 
dShen etalJl200llGuo et aT]|2009b, consistent with previous 
results (see. ICimatti et alj|2008t ICassata et al.H201d and ref- 
erences therein). 

In contrast, the quiescent galaxies associated with the z = 
1.62 cluster show a relative lack of compact galaxies com- 
pared to the quiescent field galaxies at z = 1.6 at fixed mass. 
Quantitatively, the quiescent cluster galaxies with masses > 
3 x 10 10 M Q and V z > 0.5 have an interquartile (25-75%- 
tile) range of R e ff = 1.2-3.3 kpc with a median of 2.0 kpc, 
whereas the field galaxies have an interquartile range, R e ff = 
0.9 -2 .4 kpc with a median of 1.3 kpc. The size of a typical 
massive, quiescent galaxy in the z = 1 .62 cluster is larger com- 
pared to field galaxies. This trend is c onsist ent q ualitatively 
with r ecent findings by ICooper et al.l d201 lb and IZirm et all 
d201 lb . We note, however, that the difference between the 
sizes of the cluster and field quiescent galaxies declines at 
higher masses, as many of the galaxies with M > 10 11 M Q in 
both the field and cluster samples have larger effective radii 



(> 2 kpc; similar to the fin dings of IRettura et al. 12010b . We 
note that recent work from iRaichoor et al.l (12012 ) concludes 
an opposite trend such that early-type galaxies in higher den- 
sity regions are smaller. However, the significance of this re- 
sult is likely a consequ ence of sample selec tion and analysis 
method, as discussed in lCooper et al.l (1201 lb . 

The relative lack of compact quiescent galaxies in the z = 
1 .62 cluster is unlikely a result of selection effects. There are 
inherent biases and systematics in the measu rement of both 
the effective sizes and stellar mas ses (see, e.g.. lPapovich et aD 
1200 1L 120061: lHaussler et"aTll2007b . However, these mostly af- 
fect comparisons between samples of galaxies taken from dif- 
ferent datasets and at different redshifts. In the case here, both 
the cluster and field galaxy samples are selected at the same 
redshift and using the identical CANDELS dataset. There- 
fore, the same systematics and biases affect both samples 
equally. As a result, the relative comparison between the 
galaxies in the cluster and field is robust. 

Figure [5] shows the distribution of effective radii for the 
z = 1.62 cluster and z = 1.6 field quiescent galaxies rela- 
tive to the l ow-re dshift relation for early-type galaxies from 
IShen et alJ (12003b . The CANDELS UDS samples include 
quiescent galaxies with stellar masses > 3 x 10 10 M Q and 
V z > 0.5. The main difference in the samples is that the 
cluster galaxies at z = 1.62 have a relative lack of quiescent 
galaxies with low effective sizes compared to the field sample, 
as discussed a bove. Formally, a Man n-Whitney-Wilcoxon 
rank-sum test ( iMann & Whitnevl [l947l) gives a 12% likeli- 
hood (~ 1.2cr) that both the CANDELS z = 1.62 cluster and 
Z = 1 .6 field samples are drawn from the same parent sample. 
The significance increases to ~ 2cr if we consider a higher 
fidelity sample of cluster and field galaxies with integrated 
photometric redshift probability distribution V z > 0.65. 

There is strong evidence that the size distribution evolves 
from z = 1-6 to z ~ 0, as inferred from other studies. We test 
this by computing a likelihood that the effective sizes of qui- 
escent field galaxies at z = 1.6 and quiescent cluster galaxies 
at z = 1.62 have the sam e mean size s as th e local sample of 
early-type galaxies from IShen et all (120031) . A Student's t- 
statistic gives likelihoods of 2 x 10~ 7 and 6 x 10~ 3 , for the 
Z = 1.6 field and z = 1.62 cluster samples respectively, where 
the higher likelihood for the cluster is a result of the fact that 
the mean size of the cluster galaxies is larger than that of the 
field galaxies. Regardless, based on these tests both the qui- 
escent galaxies in the cluster and field at z = 1.6 show strong 
size evolution compared to the field. 

4.2. The Evolution of the Size-Mass Relation in Clusters 

Figure [6] shows the evolution of the relative mean sizes for 
quiescent galaxies for CANDELS in the z = 1.62 cluster and 
in the z ~ 1.6 field from the compared to other samples in 
the literature. In all cases we me asure the size rela tive to the 
low-redshift relation from SDSS IShen et al.l (12003b . The fig- 
ure includes the mean relative sizes of early-type galaxies in 
other high-redshift cluster samples, including MS 1054-03 at 
z = 0.8 3 an d RXJ0152.7- 1 357at z = 0.83 from lBlakeslee et all 
J2006h and lHolden eFall (I2009I). a nd RDCS 1252.9-2927 at 
z= 1.24 from IRettura et alJ (120101) . In addition, the figure 
shows the mean relative sizes for the field sample of early- 
type galaxies of ICimatti etaT] (12008b at z ~ 1.4-2.0. The 
shaded curve shows th e best-fit relation t o the evolution of 
quiescent galaxies from lvan derWel et alJ d2008h . 

We parametrize the size evolution in figure|6]for the cluster 
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Figure 4. The left panel shows the relation between the circularized effective radii and stellar mass for quiescent galaxies in the z = 1 .62 cluster for galaxies with 
projected distances S pro j < 1.5 Mpc from the cluster center and V- > 0.5. The right panel shows same relation for z = 1.6 quiescent galaxies in the field, selected 
in the same way as the cluster galaxies but with R„ m j > 3 Mpc. In both panels the size of the data point (boxes) scales with, as indicated in the legend of 
the left panel. The unfilled boxes deno te objects detecte d at 24 /im with f v (24fim) > 50 fily. The solid and dotted lines show the z = 0.1 size -mass relation 
for ear ly-type galaxies from the SDSS I Shen et al. 2003). In each panel the filled stars correspond to the 1 < z < 2 early-type galaxy sample from Cimatti et al] 
I2008B . Quiescent galaxies in the field at z ~ 1.6 in the CANDELS data have sizes similar to these other studies. There is a relative lack of compact quiescent 
galaxies in the cluster compared to galaxies in the field. 

the very compact, pa ssive galaxies at 1 < z < 3 (see also, 
iDamjanov et al.ll201 ll) . The data for the quiescent galaxies 
in the z = 1.62 cluster suggest that quiescent galaxies in the 
high density region of clusters experience slower size evolu- 
tion from z= 1 -6 to z ~ compared to the field. 



samples as j? e ff oc (1 +z)°. Fitt ing t he data points for the clus- 
ters of (IB lakeslee etal .1120061) and iRettura et al.l ( 120101) , and 
the z = 1 .62 cluster from the CANDELS data, we find a = 
-0.6 ±0.1. This is highly consistent with the 30% increase 
in the si zes of brightest clu ster galaxies at z = 1 to 0.25 mea- 
sured b y lStott et all (1201 11) . In comparison. Ivan der Wei et"aTI 
(120081) derived a steeper exponent, a = —1.0 ±0.1 consider- 
ing samples of field and cluster early-type galaxies, and this 
rapid evolution seems required to match the mean sizes of 
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Figure 5. Comparison of the distribution of circularized effective radii of 
quiescent galaxies in the z = 1.62 cluster and the z = 1.6 field in the CAN- 
DELS UDS dat a. The size distri bution is measured relative to local early- 
type galaxies of Shen et al. 12003)). The v ertical dotted lines show the scatter 
about the mean relation (dashed line) fromfShen et an <2003h . The histograms 
show the distribution for both the z= 1 .62 and cluster and z = 1 -6 field quies- 
cent galaxies with solar masses > 3 X 10 10 Mq and V z > 0.5. The average 
relative size of a quiescent galaxy is larger compared to quiescent field galax- 
ies. This is primarily due to the lower number of compact quiescent galaxies 
in the cluster. 



5. EVIDENCE FOR STELLAR DISKS IN QUIESCENT CLUSTER 
GALAXIES 

5.1. Ellipticity Distributions 

As discussed in § 3, the quiescent galaxies in the z = 1.62 
cluster have concentrated, spheroid-dominated morphologies 
(Sersic indices n > 2). In addition, many of the quiescent 
galaxies have low axial ratios, corresponding to high elliptic- 
ities (e = l—b/a where b/a is the ratio of the semi-minor to 
semi-major axes from GALFIT, see § 2.3). Indeed, many of 
the quiescent galaxies in figure |2] show elongated morpholo- 
gi es with significant ellipti city. 

Ivan der Wei et alj (1201 ll) recently cited the high ellipticities 
of a majority of field quiescent galaxies at z ~ 2 as evidence 
that these galaxies have prominent disk components, consis- 
tent w ith other studies dMcGrath et al.l l2008t IWeinzirl et al] 
1201 ll) . This is similar to the observation for both the cluster 
and field galaxies at z = 1.6 here. 

Figure Q shows the measured ellipticities for the quiescent 
galaxies in the z = 1 .6 field and the z = 1 .62 cluster as a func- 
tion of stellar mass. Both the quiescent galaxies in the z = 1 .62 
cluster and z = 1 .6 field in the CANDELS data have rela- 
tively high measured ellipticities. The median ellipticity for 
both samples e me d = 0.4. As illustrated in figure |7l there is 
no strong evidence that the ellipticity distributions differ be- 
tween the z = 1 .6 field and z = 1 .62 cluster samples. A Mann- 
Whitney-Wilcoxon rank-sum test finds no statistically signifi- 
cant difference between the ellipticity distributions for the two 
samples: we are unable to reject the hypothesis that they are 
drawn from the same parent sample. 

The measured ellipticities of the quiescent galaxies at z = 
1.6 are comparable to the ellipticities measured for lenticu- 
lars and early-type spirals in clusters at < z < 1, which have 
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Figure 6. Comparison of the evolution of the effective radii for early-type 
and quiescent galaxies. The large squares show the relative sizes of CAN- 
DELS quiescent galaxies in the z = 1 .62 cluster (solid black datum) and the 
Z = 1.6 field (solid gray datum). The ellipses show the mean relativ e sizes 
of early-type galaxies from two clusters from Blakeslee et al. (2006, B106) 
using the stellar mass rela tion from Holden et al. (20Q2) at z = 0.83, and 
from one cluster at z = 1.24 Rettura et al. 12010, RelO). The sta rs show pas- 
sive galaxies in the field at 1.2 < z < 2.0 ICimatti et alj|2008l Ci08). All 
data points show the galaxies sizes relative to the local size-mass relation of 
Shen et al. 12003 ). The horizontal lines show this local relation its scatter. 
T he shaded curve shows th e size evolution in early-type galaxies measured 
bv lvan der W ei et al. 1 200|) for a mix of field and cluster galaxies. The thick, 
solid line shows the size evolution measured here for the cluster galaxies only. 
The fit sugg ests milder size evolution from < z < 1.6 for cluster galaxies 
compared to[yan der Wei et al. 1 20081) . 

e m°d = 0-4-0.5(e.g„ IVuicani et alJUoTlh . In contrast, mas- 
sive cluster e llipticals at low redshift h ave lower ellipticities, 
e ~ 0.2-0.3 (IHolden et al.ll2009il201 11) with no indications of 
evolution (e.g.. IVuicani et al.H2.Ql 11) . 

However, unlike galaxy samples at lowe r redshifts (e.g., 
iHolden et al.ll2Q09b Ivan der Wei et alj|2009b|) , we find no ev- 
idence for a trend between the ellipticity and stellar mass in 
the z — 1.6 field and cluster samples. As illustrated in fig- 
ure [7] the ellipticities of the higher mass (> 10 n M Q ) galaxies 
in both the z = 1 .6 field and z = 1 .62 cluster remains high, with 
a median e me( j = 0.4 and with an interquartile rang e spanning 
e = 0.3-0.7. In contrast. Ivan der Wei etafl d2009bl) find a me- 
dian ellipticity of e ~ 0.2-0.3 for non-star-forming galaxies 
at 0.04 < z < 0.08 with M > 10 11 with no apparent evo- 
lution to z ~ 0.6-0.8 ( IHolden et all 1201 lh . Therefore, the 
massive (> 10 11 Mq) quiescent galaxies in both the z= 1.62 
cluster and z = 1 .6 field have higher ellipticities than lower 
redshift (z < 1) counterparts. 

5.2. Surface Brightness Profiles of Quiescent Cluster 
Galaxies 

These ellipticities may indicate flattened disk-like struc- 
tures viewed in projection. Roughly 50% of the cluster sam- 
ple, and 30% of the field sample have e > 0.5 (see figure [7}. 
Assuming inclination angles are distributed randomly, this 
implies that a large portio n of the massive qu iescent galax- 
ies have disk components (Lamb as et alj|1992l) . We investi- 
gate the presence of disk components by studying the surface 
brightness profiles of the four most massive quiescent galax- 
ies associated with the z = 1.62 cluster, IDs 39716, 40170, 



1.0 



f 0.8 h 



II 0.6 

UJ 

0.4 
} 0.2 h 

UJ 

0.0 



1 — i — i — i — i — i — i — | 1 1 — i — 

Quiescent Galaxies with 
M > 3 x 1 10 M and P z > 0.5 

□ ■ 




□ 



HO. 2^ 
0.4 o" 



_l I d I I L_l_ 



_J I L 



0.0 



o 
on 



0.6 

Ho. 8 

1.0 



X 

< 



0.2 0.5 1 2 5 
Stellar Mass [10 11 M Q ] 

Figure 7. The distribution of galaxy ellipticity, e = (1— b/a), as a function 
of stellar mass for the quiescent galaxies in the z = 1.62 cluster (filled boxes) 
and in the z = 1.6 field (lightly shaded boxes) in the CANDELS UDS data. 
There is no statistically significant difference between the distributions for 
the cluster and field samples. The lines show the median and interquartile 
(25-75%-tile) values in bins of 0.25 dex in mass for the combined field and 
cluster sample. The median ellipticity, e ~ 0.4, at stellar masses < 10" Mq 
is si milar to values found for no n-starforming galaxies in SDSS at 0.04 < z < 
0.08 I van der Wei et al. 2009b ). However, there is no strong trend between 
the ellipticity and stel lar mass, which contrasts with observations of lower 
redshift galaxies (e.g. JHolden et al.l2009tlvan der Wei et alJ2009bt) . 

40640, and 42952 (see Table 1). These four galaxies all have 
stellar masses > 1 x 10 11 Mq (see figure [H with low levels 
of star formation^ Three of these four most massive galaxies 
(39716, 40640, 42952) have e > 0.4. 

Figure [8] shows the one-dimensional surface brightness 
profiles for these galaxies. We fit each galaxy using three 
models. These include a best-fit GALFIT model using a sin- 
gle component with the Sersic index, n, as a free parameter. 
We also considered a model with a single component with the 
Sersic index fixed at n = 4. Lastly, we considered a model 
with two components, where the Sersic index is fixed at n = 4 
for one component and at n = 1 for the other component. In 
addition, objects 39716 and 40170 show indications of faint 
companions with angular separations of less than one arcsec- 
ond. For the analysis here, we masked the light from these 
objects to prevent them from affecting these surface bright- 
ness fits. However, we find that masking these faint objects 
changes the derived effective sizes and ellipticities by <15%. 

In all cases the single component fits require Sersic indices 
n > 2 for these objects. They are spheroid dominated. Gener- 
ally, the two component models have lower residuals between 
the model and the data, particularly at larger radii (see the bot- 
tom panels for each galaxy in figure[8]l. The disk exponential 
scale length for the n = 1 components range from 2-5 kpc, 
consistent with the disk scale lengths for low redshift g alaxies 
of comparable stellar mass in SDSS dFathi et al.ll2010h . In all 
cases the n = 1 components have ellipticities that are in within 

3 We find from the analysis of their spectral energy distributions (§ 2.2 and 
3) limits on the SFRs of < 10 Mq yr" 1 , with the exception of 42952 which 
is consistent with \t < 40 Mq yr~' . These SFRs are consistent with the limits 
from their (lack of) detected Spitzer 24 (im emission, /„(24/im) < 40 /iJy, 
implying \E'(24/im) < 5 Mq yr~' . Including both the constraints from the 
24/im data and analysis of spectral energy distribution, the specific SFRs for 
these galaxies are very low, < 5 X 10~ 2 Gyr~' . 
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Figure 8. Surface-brightness profiles for the four most massive, quiescent galaxies (M > 1 X 10" Mq ) associated with the cluster at z = 1 -62, with ID numbers 
as labeled. The box points in the top panels show the measured surface-brightness profile. The inset images show a 6" X 6" cutout of the galaxy using the WFC3 
F125 and F160 data. The curves in the panels show different model fits to the F125W data. For each galaxy, the left panel shows a model with two-components, 
where one model has a fixed Sersic index, n = 4, and the other has a fixed Sersic index, n = 1 . The middle panel has a single model where the Sersic index is a 
free parameter. The right panel has a single model with fixed Sersic index n = 4. The bottom panels for each galaxy show the difference between the measured 
surface-brightness and each model. Error bars on the surface-brightness measurements are shown in the bottom panels only for clarity. 



20% of the ellipticity from the fits for each object with single- 
components. The implication is that the light profiles of these 
massive z = 1 .62 cluster galaxies are inconsistent with a model 
of constant Sersic index, and they instead favor a model with 
a radially dependent Sersic index to describe their structure, 
with high Sersic indices at small radii changing to lower Ser- 
sic indices at larger radii. 

We use the two-component models to estimate crudely 
the "bulge-to-total" (B/T) ratio for each galaxy, defined 
as the ratio of the flux in the « = 4 component to the to- 
tal flux. Objects 39716, 40170, and 40640 are bulge- 
dominated, with B/T = 0.5-0.7, consistent with those of 
lenticular galaxies and elliptica ls ( ISimien & de Vaucouleursl 
[19861 IGraham & Worievl 120081) . Object 42592 has a lower 
ratio, B/T = 0.3. The disk component dominates the 
light in this galaxy, similar to middle-to-early-type spiral s 
(ISimien & de Vaucouleursl 119861; IGraham & Worievl 120081) . 
Assuming the stellar mass traces the J\25 light, the derived 
B/T imply that 30-50% of the stellar mass lies in the disk 
component for most of these galaxies, although the disk in 
42592 may contain as much of 70% of the stellar mass. Inter- 
estingly, unlike the other three objects (39716, 40170, 40640) 
with all reside within 0.5 Mpc of the luster, object 42952 
sits at a projected distance 1.2 Mpc, and there is some evi- 



dence quiescent galaxies at this distance exhibit more disk- 
dominated morphology (R. Bassett et al. 2012, in prep). 

We conclude that these massive quiescent galaxies asso- 
ciated with the cluster show evidence for prominent stellar 
disk components in their surface brightness profiles. The 
lack of significant star formation in these galaxies suggests 
that these disks are primarily stellar systems. However, it 
may also be the case that the spheroids form through the mi- 
gration of stars formed from violent instabilities in the disk, 
which stabilizes the gas in the disk agai nst further instabili- 
ties that would oth erwise form stars (e.g. jDekel et alj|2009at 
iMartig et ail 120091) . Currently, the only i/Sr-quality data 
(FWHM ~ 0.1-0.2") for these galaxies is the CANDELS 
F125W and F160W imaging used here, and we are unable to 
test for surface-brightness gradients indicative of variations 
in the stellar populations of these possible bulge and disk 
components. However, we see no measurable color gradients 
in WFC3 J125 -H160 images. While this is consistent with 
iMcGrath et all d2008l) who find negligible ACS / 8 i 4 - NIC- 
MOS H\ 6Q color gra dients of early-type g alaxies at z ~ 1.5, 
iGuo etalJ J2011I) and lSzomoru et al.l ( 120111) find evidence that 
some quiescent galaxies at z ~ 2 have negative color gradi- 
ents with bluer cores and redder outer regions. To test for 
color-gradients in the galaxies here will require data into the 
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rest-frame near-UV using, e.g., HST/ACS observations. 

If this result generalizes to the full sample of quiescent clus- 
ter galaxies, the ellipticity distribution in figure [7] provides 
evidence that a large fraction of these galaxies host stellar 
disks. The lack of color gradients in WFC3 /i25 _ ^i60 im- 
ages suggests the stellar populations in the di sk and spheroid 
components are fairly homo geneous (see, e.g JPapovich et ail 
120051; iMcGrath etafl 120081) . However, given that massive 
cluster galaxies at low redshift are dominate d by spheroids 
with n > 4 and show no evidence for disks dPostman et al.l 
12005b iHolden et al.ll2009f) . these disk structures must be de- 
stroyed at some later time. We discuss the implications of this 
evidence below. 

6. DISCUSSION 

6.1. Accelerated Evolution in High Density Regions 

The quiescent galaxies in the z = 1 .62 cluster and the z = 1 .6 
field share many common properties. Their morphologies 
show dominant spheroidal components. However, both the 
field and cluster samples have broad ellipticity distributions 
(figure|7]and § [5j, suggesting the presence of disks. Based on 
the modeling of the galaxies spectral energy distributions and 
(lack of) IR emission, these quiescent galaxies have low levels 
of star formation (see § 2.2 and 3), implying they either have 
low cold-gas fractions, or that the cold gas in the galaxies is 
stable against i nstabilities, perhaps as a result of the domi- 
nant spheroids (iDekel et al.ll2009atlMartig et al. || 2009l). or the 



domi nance of a stellar component in the disk (ICacciato et al.l 
l20Tlh . 

The main difference between the cluster and field quiescent 
galaxies at z = 1 .6 is the relative lack of compact, massive 
quiescent galaxies in this cluster compared to those in t he fi eld 
(figure|Ui, and this is significant at the > 90% level (§ 14. U . If 
correct, then this result implies that the quiescent galaxies in 
the z = 1 -62 cluster have experienced accelerated size growth 
relative to the quiescent galaxies in the field at z = 1 .6. 

One possibility is that prior to z = 1.6 the quiescent galaxies 
in the cluster experience accelerated spheroid growth associ- 
ated with the migration of stellar clumps originally formed 
in the galaxy disks. Theoretical considerations and cosmo- 
logical hydrodynamical simulations predict that galaxies at 
z > 2 form stars from gravitational instabilities in the disks 
fed by col d gas accreted in str e ams and minor merge rs along 
filaments (IDekel et al.l l2009bt ICeverino et ail 120101) . Over 
~ 0.5 Gyr the instabilities and clumps in the disk migrate in- 
wards and merge and form a passive spheroid (IDekel et al.l 
l2009at iBournaud et al.l 1201 ll) . The spheroid stabilizes the 
disk (leaving it intact) against instabilities, which suppresses 
star formation. One feature of this model is that spheroid- 
dominated galaxies at z ~ 1.5-2 should show disk compo- 
nents with a scale length comparable to that of the spheroid 
dMartig et al.H2.Q09h . This is qualitatively consistent with our 
CANDELS observations of quiescent galaxies in both the 
z = 1 .62 cluster and the z = 1 .6 field. Nevertheless, because 
the sizes of the quiescent galaxies in the z= 1 .62 cluster are 
larger on average compared to those in the z = 1 .6 field, this 
would imply that processes of star-formation in the disk and 
the migration of stars into a central spheroid happen earlier 
and/or occur at an accelerated rate for the cluster galaxies, 
possibly as a result of enhanced gas accretion associated with 
the higher density region. 

Parenthetically, these simulations predict that a small frac- 
tion of the gas that flows into the spheroids fuels the growth 



of supermassive black holes (SMBH). IBournaud et al.l (1201 ll) 
predict that at z ~ 2 a galaxy with 10 11 M Q will have a SMBH 
with an accretion rate that corresponds to an X-ray luminosity 
of 10 42-43 erg s" 1 . One of the most massive quiescent galax- 
ies in the z = 1.62 cluster, ID 39716 (see Tableland § 5.2), 
has a stellar mass 1.6 x 10 11 M Q and has an X-ray luminos- 
ity ~ 3 x 10 42 erg s" 1 dPierre et al.ll201 ll) . Neither the near-IR 
spectrum (lTanakaetal.1120101) nor IRAC colors of this galaxy 
show any indication of an AGN. Nevertheless, if this X-ray 
flux stems from accretion onto a SMBH, then it is consistent 
with the model predictions. 

There are also effects associated with the higher density 
region of the cluster that could influence the galaxies' mor- 
phological evolution. Interactions between the galaxies and 
the intracluster m edium (ICM) play a significant role in mas - 



play . 

sive clusters (e.g., lMcCarthv et al.ll2008TlBalogh et alj|2009l) 
However, it seems doubtful that these influence galaxy evo- 
lution in this z = 1.62 cluster, as observations of its X-ray 
emission show tha t the hot ICM gas has not developed fully 
dPierre et al.l 1201 ll) . Therefore, effects associated with in- 
teractions with the ICM are likely less important drivers of 
galaxy evolution i n this cluster (see further discussion in 
IMcGee et al. l l2009l) . 

Galaxies associated with the cluster likely experience an ac- 
celerated merger rate for the reason that this cluster is still 
forming and has a high density of galaxies. Galaxy assem- 
bly via mergers is most effective in small groups and forming 
clusters at lower redshifts, because these sys t ems have lower 
veloc i ty dispersions (see , jTran et alj 120081: Mcintosh et alj 
l2008t IMcGee et al l 12009]: IWilman et al.l 120091) . For exam- 
ple. lMcIntosh et a" ( 2008 ) find evidence that near-equal mass 
("major") mergers between red galaxies are more common in 
z < 1 groups than in massive clusters. It follows that mergers 
are also an important assembly mechanism in forming clus- 
ters at higher redshift, and that this process is accelerated in 
the higher density regions. 

Mergers are expected to be an important assembly mecha- 
nism for massive, morphologically early-type galaxies. Mod- 
els show that at late times (z < 2) these galaxies grow pri- 
marily through dissipationless minor mergers and through the 
steady accretion of s maller stellar systems formed outside the 
galaxies' vi rial ra dii (lLoeb & Peeblesl200llCiotti et al.l2007t 
iNaab et al.ll2009t lOser et al.ll2010h . These events increase the 
galaxies' effective radii with a relatively mild increase in stel- 
lar mass. Measurements of the galaxy merger rate in high den- 
sity regions (such as clusters) show that these me rgers occur 
mainly without star formation (Ellison et al. 2010). (Mergers 
involving even small amounts of star formation are disfavored 
by the measured evolu tion of the colors of cluster galaxie s 
down to lower redshift. Ivan Dokkum & van der Marelll2007l) . 
A higher incidence of these dry mergers is expected to play a 
dominant role in the e volution of quiescent gal axies at z < 2 
(see, e.g., discussion in Ivan der Wei et al.ll20ll . 

Our observation that quiescent galaxies in the z = 1 .62 clus- 
ter have larger sizes could be related to an accelerated dry 
merger rate associated with the forming cluster. There is some 
evidence to support this as the massive quiescent galaxies in 
the z = 1 -62 cluster appear to show a higher frequency of com- 
panions than those in the field, which implies a higher current 
rate of mass growth from mer ging (lLotz et al.l 1201 ll) . This 
is consistent with the models of Shank ar et al.l (1201 ll) . which 
predict that at fixed stellar mass central galaxies in larger mass 
halos have large effective sizes compared to central galaxies 
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in lower-mass haloes. In the Shankar et al. (2012) models, 
galaxies in different halos masses undergo different types and 
numbers of mergers, consistent with the results here. 

While the observations suggest that cluster galaxies have 
experienced an accelerated history at redshifts greater than 
z = 1.6, there is evidence that additional evolution is also 
required. First, the most massive galaxies in the clus- 
ter are still only ~ 10-50% as m assive as the bright- 
est galaxies in low redshift clusters fe lakeslee et al.l 120061: 
iHolden et al.l 12009b IValentinuzzi et al.ll2010h .~ These galax- 
ies need to increase both their stellar masses and their ef- 
fective sizes by at least a factor of 2. Simulations pre- 
dict that the growth of massive cluster galaxies at "late" 
times (z < 1.5) should occur more through the dissipa- 
tionless mergers o f relati vely massive progenitors (e.g ., 
De Lucia & Blaizotl 12007b IRuszkowski & Springel 120091) . 
Ruszkowski & Springell pool " predict that dry major merg- 
ers are an important growth mechanism for galaxies at z < 1.5, 
and that the number of major mergers declines strongly with 
galaxy mass. If this is the case, then we may expect that the 
more massive quiescent galaxies (> 10 11 M@) associated with 
this cluster will experience ~1— 4 additional major dissipa- 
ti onless mergers (see also § 6 .2). This is similar to the findings 
of Ivan Dokkum et al.l JT999h that 50% of massive early-type 
massive galaxies in the z = 0.83 cluster MS 1054+03 will un- 
dergo a dry major merger at z < 1 . 

Second, the ellipticity distributions of the quiescent galax- 
ies in the z = 1.62 cluster are shifted to relative high values, 
and these galaxies show evidence for extended disks. Both 
facts contrast strongly with observations of early-type clus- 
ter galaxies at z < 1 (see § 5). Mergers would account both 
for the required evolution in size, ellipticity, and mass as they 
transform the sur face brightness profiles toward higher Sersic 
indices (see, e.g.. lNavarrolll990l) . 

Lastly, measurements of the luminosity function of red- 
sequence cluster galaxies show that the bright (massive) end 
is con sistent with passive evolution since z < 1 ( iRudnick et al.l 
120091) . However, a preliminary analysis shows that in the 
z = 1.62 cluster here the bright end of the red-sequence- 
galaxy luminosity function is not fully formed, and the mas- 
sive galaxies require additional mass growth mostly through 
dry m ergers without substantial star formation IRudnick et alj 

(HoH. 

Therefore, we conclude that the quiescent galaxies asso- 
ciated with the z = 1.62 cluster require additional growth 
through dry mergers to match the properties of early- 
type massive galaxies in lower redshift cluste r s. Th is is 
consistent with the find ings of iMcGrath et al.l J2008I) and 
Ivan der Wei et al.l (1201 ll) for quiescent field galaxies at z ~ 2. 
However, as discussed in § 4.2, the quiescent galaxies in the 
Z = 1 .62 cluster appear to require less size growth from z ~ 1 .6 
to z ~ compared to field galaxies in order to them to be con- 
sistent with the siz e-mass relation for quiescent galaxies in the 
field and clusters (IWeinmann et al.l r2009t IValentinuzzi et al.l 
120101) . Because dissipationless mergers of low-mass com- 
panions ("minor" mergers) produce more size growth relative 
to stellar-mass growth. If quiescent galaxies grow primarily 
through this mechanism, then it seems to follow that the quies- 
cent galaxies in the z = 1 -62 cluster will experience additional 
mergers weighted toward more massive progenitors (i.e., they 
will experience more major mergers), compared to quiescent 
galaxies in the lower density field. 



6.2. The formation of the brightest cluster galaxy? 

One interesting possibility is that the most massive qui- 
escent galaxies in the z = 1.62 cluster could merge into 
the brightest cluster galaxy (BCG). IVaientinuzzi et alj d2010t) 
show that while early-type galaxies in local clusters fol- 
low the size- mass relation of other ( field) early-type galax- 
ies (see also IWeinmann et al.l 120091) . BCGs are often out- 
liers, ha ving significantly larger effective radii for their stel - 
lar mass dBernardi et al.l20 07: Ruszkowski & Springell2009h . 
although see lLauer et al.l (120071) and lGuo et all (120091) for al- 
ternative interpretations. The most massive galaxies in the 
z= 1.62 cluster have stellar masses w2x 10 11 M Q , and these 
would require at least a factor of 2 growth (and as much as a 
factor of 5) to achieve the st ellar mass of the BCGs in lower 
redshift clusters measured by IValentinuzzi et"aT] (120101) . 

The two most massive galaxies in the z = 1.62 cluster, ID 
39716 and 40170 (see figure |8), are both near the core of 
the cluster (each within a physical distance of <70 kpc), and 
they have a projected physical separation of 126 kpc. As- 
suming these galaxies have relative velocities of > 100 km 
s -1 (about one third the estimated velocity dispersion), they 
would experience a first pass encounter in < 1 Gyr. It there- 
fore seems entirely likely that these galaxies will merge by 
z ~ 1.2. This is consistent with simulations that predict the 
progenitor of the BCG should exp erience 1-2 major mergers 
between z ~ 1.5 and z ~ 1 ( e.g.. Ibe Lucia & Blaizotl 12007b 
IRuszkowski & Sprin gel 2009). These galaxies currently have 
stellar masses 1.5 x 10 11 M and 2.1 x 10 11 M Q , and effective 
radii, 2.6 and 3.2 kpc, respectively. Assuming they will merge 
with no additional star formation, the remnant will have a stel- 
lar mass M > 3 x 10 11 M , with a more compact morphology 
(Sersic index n ~ 4), and grow in effective radiu s to >6 kpc, 
based on arguments from the virial theorem (see lNipoti et alJ 
120031) . Additional growth through mergers and accretion (in- 
cluding accretion through dynamical friction of other galaxies 
in the cluster potential) would increase the size and mass fur- 
ther, shifting the new galaxy along (or even above) the size- 
mass relationship in figure [4] (consistent with some low red- 
shift BCGs, see Valentinuzzi et al. 2010). Therefore, it seems 
we are witnessing the progenitors of the BCG in this cluster 
before they merge. 

7. SUMMARY 

In this paper we discussed morphological properties of 
galaxies in a z = 1.62 cluster XMM-LSS J02 182-05 102 using 
partial near-IR coverage from //5T/WFC3 as part of CAN- 
DELS. The cluster shows a prominent red sequence domi- 
nated b y galaxies with colo rs consistent with passive evo- 
lution dPapovich et al.l 120101) . although there is a popula- 
tion of star-form ing galaxies in this cluster with high SFRs 
dTran et alJl2010h . Recent Chandra data for this cluster show 
that the X-ray emission is mostly attributed to point sources, 
suggesting that thi s cluster is still in the process of collapse 
dPierre et alj|201 ll) . Therefore, we are able to study galaxy 
evolution in the high density region of a forming cluster at 
high redshift. 

The HST/WFC3 images show that the cluster galaxies 
exhibit a clear color-morphology relation, where galaxies 
with colors of quiescent stellar populations have dominant 
spheroids, and galaxies with colors consistent with ongoing 
star formation have disk-like and irregular morphologies. 

The quiescent cluster galaxies follow a size-mass relation- 
ship, but the cluster is deficient in quiescent galaxies with 
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compact effective radii compared to quiescent galaxies in the 
field at z = 1.6. The average effective radii of the quiescent 
galaxies in the cluster are larger compared to quiescent galax- 
ies in the field at fixed stellar mass (> 90% significance). 

If the difference in effective radii between the cluster and 
field galaxies is generalizable, then it implies that the qui- 
escent cluster galaxies experience accelerated size evolution 
at redshifts greater than 1.6 compared to similarly selected 
field galaxies. Because other mechanisms associated with in- 
teractions between the galaxies and the cluster 1CM are not 
yet operating, we argue that to explain the observations quies- 
cent cluster galaxies have had accelerated spheroid formation, 
possibly as a result of the migration of stars formed in disks, 
and/or merger histories associated with the formation of this 
cluster. This gives rise to an accelerated size growth compared 
to galaxies in the field. 

The morphologies of quiescent galaxies in the field and 
cluster are dominated by spheroids. However, their elliptic - 
ity distributions are broad, with median values e mec j = 0.4, 
with no trend between ellipticity and mass, in contrast to 
lower redshift samples. Both the ellipticity distributions and 
the surface-brightness profiles of the massive cluster galaxies 
suggest these galaxies host stellar disk components. Because 
the quiescent galaxies have low estimated SFRs, these disks 
are likely composed of stellar material with low gas fractions 
available for star formation, either because they have depleted 
their gas supplies, or that the dominant spheroids stabilize 
the gas in the disks, hindering the formation of instabilities. 
This is true even for the massive quiescent cluster galaxies 
(M > 1 x 10 11 Mq), which show no indications of star for- 
mation, are spheroid-dominated, yet show extended emission 
consistent with disks of scale lengths, 2-5 kpc. The implica- 
tion is that these galaxies have significant stellar disks, similar 
to the interpre tation of data for pa s sive galaxies in the fiel d 
atz~ 1.5-2 (iMcGrath et al.ll200l Ivan der Wei et al.ll201 ll) . 
These extended disks are not pre sent in quiescent ga laxies in 
clusters at lower redshifts (z < 1) dHolden et al.l2009l) . There- 
fore, some mechanism must transform or destroy these disks 
in the few billion years from z ~ 1 .6 to z ~ 1 ■ 

The quiescent galaxies in the cluster at z = 1 .62 require ad- 
ditional growth to match the observed properties of massive 
galaxies in clusters at lower redshift. Several lines of evidence 
suggest this additional growth occurs via dissipationless (dry) 
mergers. These mergers will increase the sizes and stellar 
masses of the quiescent galaxies, and affect the morphological 
transformation to more compact surface-brightness profiles 
(n ~ 4). However, because the quiescent galaxies associated 
with the cluster at z = 1 .62 have larger sizes, they appear to re- 
quire slower size growth at later times (z < 1 .6) compared to 
galaxies in the field. The evidence for this is a result of com- 
paring the data here with results from the literature. The size 
evolution of massive cluster galaxies is relatively slow from 
z ~ 1.6 to the present, with sizes evolving as (1 + z )-°- 6 ± 01 
comp ared to « (1 + z) _1 for field galaxies (Ivan der Wei et al.l 
|200i. 

To summarize, the data provide evidence that the morphol- 
ogy and size evolution in the quiescent cluster galaxies is ac- 
celerated compared to field galaxies prior to z = 1.6 to ac- 
count for the larger average sizes of the quiescent cluster 
galaxies at this redshift. In addition, we conclude that ad- 
ditional growth is necessary for these galaxies and that most 
of the growth occurs via dry mergers without significant star 
formation. Furthermore, in the case of the quiescent cluster 



galaxies at z = 1 .62, much of this merger activity must oc- 
cur between 1 < z < 1.6 such that these galaxies have the 
morphological properties (Sersic indices and ellipticities) of 
cluster galaxies at lower redshift. This merger scenario ap- 
pears consistent with semianalytic model predictions, which 
predict that dissipationless mergers dominate the mass growth 
of massive galaxies at z < 1 .5 (e.g., lbe Lucia & Blaizotll2007t 
iRuszkowski & Springef2009l) . 

One caveat to these conclusions is that our analysis is 
based on only a single cluster at z = 1 .62, which currently 
has H5T/WFC3 imaging for ~ 20 quiescent cluster galaxies. 
Clearly, extending this analysis to galaxies in other clusters at 
z > 1 .5 is required to determine if the results here are general- 
ized to quiescent galaxies in other high density environments. 
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Table 1 

Properties of z = 1.6 Galaxy Samples in the UDS CANDELS field 
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Note. — A portion of the table is shown here. (1) Object ID, (2) right ascension (J2000), (3) declination (J2000), (4) photometric redshift, (5) integrated 
photometric redshift probability distribution function, see § 2.2, (6) Suprime z, UKIDSS / color, (7) UKIDSS J, SpitzerflRAC 3.6 /im color, (8) WFC3 F125 
magnitude measured from GALFIT, (9) circularized effective radius, (10) Sersic index, (11) ellipticity, e = i-b/a, (12) stellar mass and (13) star formation 
rate from analysis of spectral energy distribution, (14) projected distance from center of the z = 1.62 cluster XMM-LSS J02182-05102. 
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